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Abstract 

^-H We explore the possibility that the dark matter is a condensate of a very light vector boson. 

Such a condensate could be produced during inflation, provided the vector mass arises via the 
Stueckelberg mechanism. We derive bounds on the kinetic mixing of the dark matter boson 
^ with the photon, and point out several potential signatures of this model. 

O 

1 Introduction 

In the past decades there has been mounting evidence that approximately 20% of the energy den- 
sity in the universe is nonbaryonic, pressureless, and very weakly interacting. Structure formation, 
baryon acoustic oscillations, galactic rotational curves, the Bullet cluster - all of these theories or 
observations point towards nonbaryonic invisible stuff, with the same equation of state as nonrela- 
tivistic matter. Unfortunately, there is no compelling experimental evidence that dark matter has 
any nongravitational interaction with the standard model. At present there are two large classes 
^ of dark matter theories which are motivated by compelling solutions to other problems beside dark 

matter: Axions and various variations on Weakly Interacting Massive Particles. Both are poten- 
tially detectable using specific techniques. It is however possible that neither is correct. Amongst 
^NJ some of the less canonical candidates lie several variants of light massive vector particles [l}|3]. 

l/-^ Quite a few authors have already considered some of the consequences of the existence of such 

O particles [3||8] . However, the conclusion of these authors is that it is difficult to obtain sufficiently 

cold dark matter from a light vector particle. 
^ We propose a variation that allows us to generate an extremely cold light vector component of 

. ^ the Universe, with a pressureless equation of state. Our inspiration comes from the fact that nature 

^ seems to make use of almost every renormalizable Lagrangian term: Abelian and Non-Abelian 

gauge theories, Yukawa couplings, (/>^ theory, and yet even though the Stueckelberg mass [9 
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renormalizable in 4 dimensions it does not appear in the standard model. If nature indeed does 
not make arbitrary choices among consistent theories then there should exist a spin one field with 
a Stueckelberg mass. Therefore, if this boson is sufficiently weakly coupled to the Standard Model, 
it is worth considering whether it is a good candidate for dark matter. It is easy to populate the 
Universe with this particle: like the axion \TV\, during inflation, the expectation value of a light 
boson fluctuates. Immediately after inflation the value of the field in our horizon is a randomly 



selected (or perhaps anthropically selected 12-141) initial condition. After infiation, when the 



Hubble constant is of order of the boson mass, the field begins to oscillate. This oscillating field 
may the thought of as a Universe-sized Bose-Einstein condensate, as described in section [3j Such a 
particle is allowed to kinetically mix with the photon via a renormalizable interaction. Therefore, at 
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some level, it presumably does mix, although no lower bound on the mixing parameter is required 
for the model to work. In section [S] we find the upper bounds on the kinetic mixing parameter such 
that the Early Universe neither thermalizes nor evaporates this condensate. We also ensure that 
the vector boson lifetime is sufficiently long, and consider constraints on the coupling from possible 
apparent Lorentz violating effects. 

2 A Model of Light Vector Dark Matter 

Our massive vector will be represented by (p^ in a Lagrangian of the form: 

_ £ = 1 ^F^-F^^ + ^/^-^^^ + 2x0^^F^.) + ^cp^^l^^ + J^A>^ (2.1) 

where and F^'^ represent the field strength of ordinary photon, is the ordinary charged 
current and (p^'^ = d^(j)^ — d^(j)^. Applying a non-unitary transformation {A ^ A — x4> and 
(/)—)• + 0{x^)) we can redefine our fields in terms of the mass eigenstates called massless photon 
and heavy photon: 

-C = \ {F^^F,, + </.'^>^,) + ^cp^r + MA" - (2.2) 

By rotating {A = A — X'P and (j) = (j) + xA) we can reach the flavor eigenstates, called interacting 
and sterile photon. These two mix through their mass term: 

-C=\ (f>^^F^, + ^z^^^,) + ^(0^ - x\){r - XA^) + Jm^^ (2-3) 

Unless otherwise stated, we will use the mass eigenstate basis. In this basis, the heavy photon 
couples to the electromagnetic current with the coupling constant scaled by e — t- x^- In the limit 
M = 0, we could perform a rotation between the degenerate heavy and light eigenstates to a new 
set of states, one of which would be massless and completely decoupled. However, if we assume 
the heavy photon is the dark matter, with a finite energy density, produced via the misalignment 
mechanism, then the number density is inversely proportional to the mass and therefore it makes 
no sense to take this limit. 

Note that the model has 2 free parameters: M and x- Fundamental theory gives us little 
guidance for their values. The theory is technically natural for any values of M and Xi in the sense 
that for a cutoff of order the Planck scale, the renormalized values are of similar size to the bare 
values. If we assume that the U{1) of the standard model is grand-unified into a semi-simple or 
simple group, then x can only be induced via loop corrections. In this case, if the mass of the 
particles in the loops mc is below the grand unification scale Aqut, the natural size of x is of 
order {g'^ / {IGt:'^))"' , where n is the number of loops required to induce the kinetic mixing, and g 
is the relevant combination of coupling constants in the loops. We will see that for cp to be viable 
dark matter, x has to be tiny, less than 10~^ over the entire mass range, so for g ^ 1, n should be 
greater than or equal to about 3. If the particles in the loops are heavier than the grand unification 
scale, there is an additional suppression of at least {Kqut / 'mc)'^ ■ 
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3 Misalignment Mechanism for Vector Dark Matter genesis dur- 
ing Inflation 

The misalignment mechanism for producing a boson condensate has been considered in connection 
with the axion [11,15-18 and various other hght scalar fields such as moduli. Spatially varying 



modes of a bosonic field will be smoothed by the expansion of the universe. However the zero- 
momentum component of the scalar field A in the FRW background has the equation of motion: 

A + 3H{t)A + m'^A = (3.1) 

which is reminiscent of harmonic oscillator, with a time dependent damping term H{t). In the 
early Universe, H(t) ^ m, the scalar is effectively massless and its Compton wavelength does not 
fit into the horizon. The field is stuck: it does not go through a single oscillation and therefore 
we observe no particles. The value of the field is assumed to take on some random nonzero value, 
because when the mass term is negligible there is no reason to prefer a field value of (p^ = 0. An 
episode of inflation will generally produce a spatially uniform field, but for m <^ H in any causally 
connected patch of the universe the mean value of the field takes on some random, non zero value. 
After inflation, the Hubble constant begins to decrease. As soon as the discriminant 9H^ — Atv? 
becomes negative, the field A begins to oscillate and we can quantize the different modes and call 
them particles. Since, up to the small perturbations in the temperature, H{r]) is everywhere the 
same, the transition happens everywhere in the Universe at the same time (in the rest frame of 
A). We are left an energy density which may be thought of as a coherent state of a macroscopic 
number of particles. The particles are extremely cold and nonrelativistic, whatever their mass. 



An adiabatic perturbation spectrum arising from the fiuctuations of the inflaton fleld 18, 19 will 
imprint adiabatic spatial variations on the density of the scalar particles, as is needed to flt the 
WMAP data. On large distance scales compared with the particle Compton wavelength 1/m the 
dynamics of gravitational structure formation is identical to that for any weakly interacting massive 
particle. 

Note that inflation will produce isocurvature perturbations arising from fluctuations of the 
scalar fleld A. Such perturbations are highly constrained, and will place an m dependent upper 



bound on the inflation scale 11 , 13 , 18 , 20 - 24 for this scenario. 

We can show that the same scenario applies to a light massive vector in a FRW Universe. As 
shown in the appendix, the equation of motion for such a vector is: 

- ir^v^) = -M^rv^ (3.2) 

As inflation blows up a small patch of space, we can assume the dark photon is uniformly distributed 
and picks a particular polarization. This means that in the Cosmic frame di(j)^ = 0, and the time 



component of (|3.2|) implies = as long as M 7^ 0. The spatial component of (|3.2|) satisfies: 



(j)' + 3H(j)i + M20* = (3.3) 

We see that each spatial component of the vector satisfies the same equation of motion as the 
scalar A in the previous example and so has the same dynamics. After entering the lightly damped 
oscillation regime the vector behaves just like dust with d{pa?)/dt = where p = {M'^cjP). Taking 
the upper bound of = nipi when M ~ we can see that the mass of (j)^ should satisfy M > 
nl,j^jHoh = 6.6 X lO^^^eV. This mass corresponds to a wavelength of about lO^^pc. This lower 
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bound on the mass is weaker than the one imphed by the existence of compact galaxies j25j with 
L ~ 1 kpc and Mcg ~ 2 x W^^Mq. Requiring that the Compton wavelength of the dark matter 
is low enough to allow structure formation on the kpc scale gives a sharper bound on the lowest 
mass: 

1 kpc < — = — ^ M > 1.67 X IQ-"^"^ eV (3.4) 

Ap MVesc ~ 

The amount of dark matter produced by this mechanism becomes simply a randomly chosen 
initial condition for the value of the field in our patch of the universe. In other regions of the 
universe, which are beyond our current horizon, the dark matter abundance is different. In ref. |14| 
it was shown, that for an axion or similar dark matter condensate produced during inflation, 
assuming other parameters do not vary, the regions of universe with dark matter abundance of order 
the abundance in our observed universe are the most highly correlated with physical features of our 
universe that seem favorable for existence of observers, allowing for an "anthropic" explanation of 
the dark matter density. 



4 Stueckelberg versus the Higgs mechanism 

The vector mass M^0'^(/>^ is not manifestly gauge invariant. In the Standard Model of particle 
physics, all massive vector particles acquire their mass due to a Higgs mechanism. However, if the 
(p were to get its mass from the Higgs mechanism, the inflationary misalignment mechanism will 
not work to produce a condensate. Assuming the Higgs Lagrangian is: 

£ = [{d^ + ig<PM^ + A((^2 _ m2/(2A))2 (4.1) 

the mass term for <p is M'^cf)'^ = g^m^cfP' / {2X), however, the symmetry breaking happens around 
~ I , which implies that the is massless above this temperature. Therefore, in order to 
make sure that there exists a time when M < H{T) while (p is not massless, we need to satisfy: 



_ 2AM 



l<H/M = ^j—_ = ^j-^_=-— (4.2) 



Therefore we need: 



M g^ , , 

> ^ (4.3 

nipi - 2A ^ ^ 

We can look at the Z boson to illustrate this condition: the right hand side is of the order 
2g^v'^/mj^ ~ 10~^ even for a heavy Higgs (500 GeV) and so a condensate of W and Z bosons 
could not have been created by a misalignment mechanism. However, one could imagine taking 
the limit in which = g'^v'^ = (p'rr? jlX is fixed, but both m/i — >• oo and A — ?■ oo. In this case the 



right hand side of (4.3) can be made arbitrarily small and the vector retains its mass for arbitrarily 
high temperature. The limit rrih — )• oo can be handled in a better way: parametrize the Higgs in 
polar coordinates = (u + /i)e*^/^ and integrate out the heavy h. The effective Lagrangian of the 
light degrees of freedom takes the form: 

£ = -^F^- ^(M^^ + 5^e)2 (4.4) 
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which is identical to the Stueckelberg Lagrangian [9, 10 , with 9 filhng the role of the Stueckel- 
berg scalar field which fixes the correct number of degrees of freedom for a massive vector. This 
Lagrangian is still invariant under: 

AxA = A + dX 

AxO = 9 -MX (4.5) 

A redefinition (f)^ = A^^ + d^9 leads to F^^^ = (f)^^ = d^^cj)^ — d'^tj)'^ and gives us a massive vector 
described by: 

Cs = -^r-cP^,-^qy^cl>^ (4.6) 

Naturally, this Lagrangian is still invariant under A;^, although it is not invariant under the naive 
gauge transformation (p^ — )• (p^ + d^X. Unlike the nonabelian case, for a C/(l) gauge theory, the 
Higgs boson is not needed to unitarize the scattering of the longitudinal mode of a massive vector 
boson, and is unnecessary for renormalizability. 



5 Bounds 

5.1 Early Universe - Compton Evaporation 

In order to be a successful dark matter candidate, the dark photon has to be a stable particle both 
in vacuum and in the dense, ionized early Universe. For light bosons, we need to ensure that 
the dark photon population does not get thermalized, otherwise it would become ultra relativistic 
and fail to be a good dark matter candidate. As with photons and plasmas, the main process for 
thermalization is the Compton- like scattering process: (j)e^ — )• (j)e^. However, this process will be 
suppressed by a factor of with respect to two other processes: (pe^ -f-)- 76^. We will call the right 
going process Compton evaporation and the left going Inverse Compton evaporation. Therefore in 
order to ensure there are enough dark photons left after interaction with plasma, we need to require 
that Compton evaporation rate F is smaller than the expansion rate of the universe H{T). Such 
condition will also imply that the thermalization rate from Compton-like scattering will be small 
and we will be left with enough cold dark matter to populate our Universe. In order to investigate 
this bound we need to know the product of the velocity and cross-section va{M,p) as a function of 
the dark photon mass M and electron three momentum p, which can be re-expressed for M <^ me- 

X 8a^Y^7r(3m2 + 2p2) 
v<t{M,p) = / ' + 0{M) (5.1) 

The width of the dark photon in plasma is then given by the thermal average over the electron 
momentum density distribution for a given temperature of the Universe. We would like this width 
to be smaller than the characteristic expansion rate of the Universe at given temperature: 



H{T) > F(r) = j dp^a{M,p)vip)n{p,T, fi{T)) (5.2) 



Where we have used the exact a{M,p), not the approximate expression (5.1), n{p,T,fi) is the 
Fermi-Dirac distribution with chemical potential We have chosen // = for T > rUe and after T 
drops below nie it was picked to be consistent with today's electron co-moving density. Given that 
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the early Universe is growing less and less dense, the strongest bound on x is in effect at the earliest 
time the dark photon is present, that is at the time when the misalignment mechanism kicks in 
at M H. This guarantees that if the dark photon survives the first characteristic time period, 
then it will not evaporate anymore during the subsequent time. The condition H ^ T does not 
guarantee this, but is a lower bound on such survival. We find it is unnecessary to consider other 
particles than the electron, since the contribution of all other charged particles with mass mj and 
charge qi will be suppressed by a factor gi{me / ruif' {qi / q^Y which together with their suppressed 
thermal momentum distributions will make their contribution small. Likewise, it is unnecessary 
to consider other evaporation processes such as cj)^ — )• 77 since they become important for dark 
photon mass of order M = {iTipim'^y/^ ~ 10^^ eV - well above the range we consider in this paper. 
The bound imposed by Compton evaporation is plotted in Figure [2] and labeled Early Universe. We 
would like to point out two features. When the Universe reaches temperatures of order T ~ O.lmg, 
its free charge density significantly drops and the evaporation process becomes much less effective. 
This temperature marks the generation of dark photon with mass M ~ T^/nipi ~ 10"^* eV, hence 
the sharp dip in the bound on x in this region. On the other hand, since the cross-section starts 
dropping off when ^/s ~ m1/M and ~ T, we can estimate a change in the slope of the bound 
around M ~ (m^/mpi)^^^ which agrees with the observed dip at M ~ 10~^ eV. 

5.2 Decays 

Apart from Compton evaporation, we can consider pure vacuum decay processes, which become 
significant once M > rue or M > Myy. Requiring that the dark photon is stable on cosmological 
timescales requires that ^ F j < i^o : 



Process 


Width 


Notes 


4> m- 

(f> —7- 777 


r, = x'a^^M^ 

3 116640007r3 \me J 


M > 2mi, Exact 

Estimate of the loop process 

See [4], valid for M < rUe 



The bounds imposed by these decays are plotted in Figure [2| 
5.3 Earth Detection 

Although (f)^ does not satisfy Maxwell's equations, its coupling to ordinary matter is the same 
as that of the photon A^. Therefore, a nonzero (/)^ will appear as a combination of electric and 
magnetic fields with strength suppressed by a factor of x- Such fields will be detectable in various 
precision experiments and it is our desire to quantify the expected phenomena as accurately as 
possible. 

By our hypothesis, in the dark matter rest frame, (j)^ = 15^3^3 cos(Aft) and so it will mimic an 
electric field E3, = xA^M sm{Mt). Given that the local density of dark matter is Tqo = M^^|/2 = 
0.3 GeV/cm^, we can infer that the amplitude of the electric field will be 

E= ^/2x 0.3 GeV/cm3/eo ~ 3300x V/m (5.3) 
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However, there is no reason to believe that the dark matter rest frame is identical with the 
Earth frame and hence we need to perform a Lorentz boost to the right frame. Given (j) = 
{(px, <Py, (t>z) cos(Mi) and v - the velocity with respect to the dark matter rest frame, the B-fields in 
the Earth frame will be: 

B = V X L^{(j)) = jMv X (fcosi'jMt) (5.4) 

We should note, that at v = 0.001c, 7 = 1 + (lO"''), and that \M(j)\ is the magnitude of electric 
field in the dark matter rest frame. Therefore the B-field is simply B = v x E - precisely as 
expected. 



5.3.1 Attenuation 

If these fields are to be detected by Earth based experiments we need to check that the dark photon 
field is not screened by the atmosphere or by the many shields that experimental physicists put up 
in order to protect their experiments from stray electric and magnetic fields. In materials, bound 
electrons will only contribute to shielding if M falls close to some energy gap of a kinematically 
allowed transition, however, even such transitions will be suppressed by factor of x^- On the other 
hand, free electrons in metals will allow a continuum of transitions, that would lead to Compton 
evaporation effects. Therefore, we will treat the interaction of dark photons with materials as a 
wall penetration by weakly interacting particles, similarly to what we have done with the early 
Universe. The change in dark photon density will be proportional to: 

a = exp ^— y dx n{x)a{M,v)^ (5-5) 

where n{x) is the free electron density of the shielding material and cr{M, Mv) is the Compton 
Evaporation cross-section and v = O.OOlc is the assumed velocity with respect to the local dark 
matter flow. Given that the respective average densities of free electrons in the ionosphere and 
copper are Uat ~ 3 x 10^^ m^^ and ncu ~ 10^^ m~^, it is clear that a whole column of 1000 km of 
atmosphere corresponds to a layer of metal about 10~^^ m thick, which is much less than any normal 
electric shielding of earthborne experiments hence we can disregard this contribution . Moreover, 
the early Universe bound on x gives an attenuation length longer than 1 meter in a copper plate 
and once combined with the bounds from the next section the attenuation length is larger than 
10^° m. 



5.3.2 Atomic Physics 

The Stark effect associated with the background dark electric field would induce a shift in the 
ground state energy of a hydrogen-like atom of order AEs = —me{2>aQeE(i/2h)'^ ~ 10~^^ eV 
which is 5 orders of magnitude smaller than the current limits [26], even if % = 1. The Zeeman 
effect would produce a shift of AEz = 5x x 10~^^ eV. This is still too small to register. The 
advantage of the Zeeman effect is that it is first order in the fields, hence in x, which makes up 
for the fact that the magnetic field is suppressed by a factor of v/c. However, effects linear in 
fields go as cos(Mt), implying a zero time average, and so a search without prior knowledge of M 
would be time consuming. However, in the region of small mass (M < 10"'^^ eV ~ 1 year) the slow 
oscillations imply no need for averaging. In this regime the slowly changing background electric 
field would mimic a slow drift in a. As an example we can take a system comprised of two clocks: 
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one driven by two photon transition from Is — )• 2s in hydrogen and the other by the hyperfine 
transition in cesium. The major correction to the hydrogen clock rate comes as a Stark effect with 
a relative shift in the frequency that goes as 

,5a; ^ A.Ei, - A^2s ^ AE2s [(2,0,0 |egz| n, 1,0) p 

UJ Els — E2s Els — E2s ~^{E2QQ — EnlQ){E2QQ — Eiqq) 

Notice that the n = 2 term dominates the sum since the degeneracy of the 2s and 2p states is 
broken by the lamb shift with /S.E{2p — 2s) ~ 10~^ eV, whereas the rest is on the order of 1 eV. 
Therefore: 

Wh {l\I^Jlamb){^ls^2s) 

Note that the Zeeman shift is identical for the Is and 2s orbitals and so there is no contribution 
linear in x- 

In cesium the Stark shift does not distinguish the states, but the Zeeman effect contributes by 
splitting the hyperfine triplet into three distinct levels, and induces a change in the clock frequency 
on the order: 

= l-5x X 10 (5.8) 

Clearly the cesium clock effect dominates for small x- Therefore, as £ oscillates very slowly, the 
experiment sees a drift in 6u}/uj which could be (naively) interpreted as drift in a of the order: 



a 



.Suj ^ _ _ « / M 



However, if the frequency of the oscillations is comparable to the time scales of an experiment, 
such as sampling rate and averaging times of individual data points, the sensitivity becomes more 
complicated. We pick 27 as a model example to illustrate our point. Fisher et al. made measure- 
ments in June 1999 and February 2003, which, given the spacing between these two dates can be 
interpretted as two measurements separated by 44 months (T = 1320 days), each averaged over 
roughly one month (to = 30 days). Therefore, the experiment should perceive a change in the value 
of the field equal to: 

S(j)(^PQ) = ^ f f dtcos[M{t + T)/h + ipo]- [ ° dtcos[Mt/h + (po]] , (5.10) 
to \Jt=o Jt=o J 

where ipo is an unknown phase of the field. Performing the integral and factorizing gives us: 

. 4<Ao . fMT\ . fMto\ . ( mto+m 

This means that for certain finetuned phases ipo ~ — M(to + T)/2 the experiment could see noth- 
ing by simply being unlucky. However, we know that 95% of time \sm{LpQ + M{to + T)/2)\ > 
sin (0.05/4 x 27r) = 0.0785 and so 95% of time 5(j) is larger than: 



> 



4(/.o . fMT\ . fMto\ . /57r\ 

sm sm sm 

Mto \ 2 J \ 2 J \200j 



(5.12) 
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We use this expression to put a 95% confidence bound on x and plot it as a-drift in Figure [2j 
Note that in the event that the sampling frequency of the expriment is a harmonic of the the 
oscillation frequency of the field, the experiment will also become insensitive to such a drift. This 
would show up as an oscillatory behavior in the bound on x and we have replaced the region where 
these oscillations become too narrow to display with a dashed line in Figure [2] 

We would like to conclude the analysis of the fine structure constant drift bounds with two 
notes. First, as the cesium contribution dominates and the exact interaction of different atomic 
levels in cesium is beyond the scope of this paper, we would like to shelve this bound as tentative 
and in need of focused treatment. Second, presence of dark matter in form of dark photon only 
mimics a drift in a and could be potentially resolved from an actual drift if one were to measure 
different energy splittings which depend on different powers of a. 



5.4 Adiabatic Conversion 



In the flavor basis, the dark and ordinary photon mix through the off-diagonal mass terms. In a 
thermal environment the mass matrix takes the form: 



M2 + 0(^2 



(5.13) 



where m^(x, t) = e^ne(x, t) /m-g is the plasma mass, which may depend on time or position. Should 
the plasma mass be slowly varying then there could be an adiabatic conversion between different 
states. Mirizzi et al. explore this effect in the context of changing electron density in the Universe 
as it expands and distorts the CMB through an excess of converted dark photons fs] , and they offer 
a very useful comparison of this process to the neutrino MSW effect. We observe that this process 
could be much more severe in the environment of ionized gas that forms a significant portion of a 
typical cluster of galaxies. 

Figure [T] shows the energy of an eigenstate of the mass matrix as a function of radial distance of 
a particle from the center of the cluster. As an example we will follow a dark photon that is infalling 
into a cluster. If the dark photon infalls adiabatically, that is slowly enough, then it stays in the 
same eigenstate of the mass matrix which in fact contains more of the original photon state after it 
crosses the point where ~ M. Therefore, the dark photon is converted into an ordinary photon, 
which thermalizes very quickly (the cluster gas temperatures are in the range 10^ — 10^ K, [28j). 
Photons generate pressure and as a result the cluster loses its gravitational glue holding it together. 
Since we do observe clusters of ionized gas, it is imperative that the section of parameter space is 
excluded. 

What does slow enough mean? In order to cross from one level to another we require that the 
characteristic time associated with the change in the system needs to be on the order of the gap 
between the energy levels. The rate of change of photon plasma mass close to the point where the 
energy gap is minimal is: 



\m^(x)=M 



1 dm-y{x) 



m^{x) dt 



V dm^{x) 



rriry (x)=M 



dx 



V n'{x) 



m-y{x)=M 



2 n(x) 



m-y{x)=M 



(5.14) 



On the other hand the mass gap between the states is minimal when M = and turns out to be: 



M=m-, 



XM 



(5.15) 
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We take a free electron density curve from 29 , replotted in[T| to determine the portion of {x, M) 
parameter space in which this infaU turns out to be adiabatic, taking the velocity of infall to be 
the escape velocity v{r) = sj GMc/r. We note that this mechanism will only work for a mass range 
of 10"^^ eV - 10"^^ eV, the lower hmit coming from the density of voids and the upper from the 
highest densities inside clusters. We have plotted the resulting region in Figure [2] and marked it 
AdC. 



5.5 Breaking Lorentz Invariance 



The existence of dark matter necessarily causes apparent Lorentz violation because it defines a 
preferred frame - its own rest frame. The effects of this frame can be measured through its coupling 
to the standard model particles. However, even if those couplings were zero, in our case this 
corresponds to x = 0, there would be a gravitational interaction. Even in the dark matter rest 
frame there is additional Lorentz violation due to the polarization of the dark photon. 

Moreover, we can see the gravitational violation of Lorentz violation by looking at the stress- 
energy tensor: assuming the polarization points in the z-direction Ai = 5i'iAcos{Mt), we get: 

/-I \ 



cos(2Mt) 



V 



cos(2Mt) 



(5.16) 



cos(2Mt) / 



The time average of Tn^ corresponds to pressureless dust, just as we concluded from 3.3 Moreover 



at late times the frequency of oscillations is shorter than the expansion rate of the Universe. The 
field begins to oscillate when M ~ ff, and at this time the oscillations cannot be averaged over a 
period and the dark matter contribution to the Stress-Energy tensor is not rotationally invariant. 
However, at that early time radiation dominates the energy density of the Universe and dark matter 
is a minor perturbation, therefore the lack of rotational symmetry of T^^, does not produce any 
significant effect. 



6 Summary 

A nonrelativistic condensate of light vector particles could be produced during inflation and is a 
viable candidate for the dark matter component of the Universe. For ultralight vector particles, 
a small kinetic mixing term with the photon could allow this particle to be detectable. After 
considering the constraints on the mixing parameter from cosmology and astrophysics, we find 
that there are some regions of parameter space which could give unusual laboratory signatures of 
dark matter, such as apparent time dependent shifts in electromagnetic properties of matter and 
dark matter conversion to visible photons in plasmas. It is peculiar that although Dark Photon 
picks a unique direction in the Universe, it lacks a mechanism to imprint this direction onto the 
Standard model contents of our Universe. This model offers a unique experimental signature - weak 
background electric and magnetic fields that cannot be screened. 
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Appendix 

A Equations of motion in the Early Universe 

The kinetic term in the Lagrangian for a massive vector L = 0'*'^ 0^1^/4 + ^1^(1)^/2, can be 
simphfied to: 



6C 



Therefore, the equation of motion in curved space reads: 

d'^iV^M = M'^cp^^ (A.l) 

Where in FRW metric this means: 

-5oa3(t)0o/3 + a\t)di(l)i,3 = a^{t)M'^(l)f} 

-a^{t) {do4>op + 3d/a(/.o/3 - a^{t)di(t>ifi + M^^^) = (A.2) 

Keeping in mind that for our candidate ^^^a = in the cosmic frame, a{t) ^ after or during 
inflation and ala = the last hne turns into: 

do{do(i)fi - d0(t>o) + 3H{do<i>p - 9/300) + M^^^ = (A.3) 
Therefore, the time component ^ = gives us: 

ao(ao0o - docpo) + SHidocpo - docpo) + M^^o = 

M^(t>o = (A.4) 

On the other hand, the spatial component P = i implies the equation for a Hubble-damped harmonic 
oscillator: 

do{do(l)i - di(f)o) + 3H{do(l>i - di(t>Q) + M''ct>i = 

dodo(f>i + mdo(t>i + M^(f>i = (A.5) 



11 




a) the s-channel b) the tt-channel 



B Compton Evaporation Matrix Elements 

For reference we have evaluated the matrix elements for the Compton Evaporation. The momenta 
were assigned as follows: 

With this convention the matrix element becomes: 

Which implies that: 

/|r|2\ - 647r2x^a^ / (m^ + M'^){Arn^ - t) - Sm^M^ 
^' 3 (m2-n)(m2-s) ^ 

m'^ + 2m^M2 + m?{3s + u) - us m'^ + 2m? M"^ + m?{'iu + s) - us\ 

{m? — s)2 irn?' — ■u)^ j 

Which agrees, up to a number of polarizations factor, with regular Compton Scattering in the limit 
M^O, x = 1. 
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1 10 100 1000 

r(kpc) 

Figure 1: Mass mixing in plasma: The solid and dashed curves show the eigenvalues of the mass 
matrix as a function of the radial position inside the cluster. The dotted Hne shows the density of 
the ionized gas in the cluster also as a function of the radial position. In order to make the level 
crossing visible we have adopted M = 10~^^ eV and x = 0-2. 
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Equivalent Scale 




Figure 2: Summary of Constraints: The early Universe behavior puts a dominant bound on x ™- 
the higher mass range, for M > 2me the bounds are dominated by decays. The Shaded region 
called AdC marks the possible combinations of (x, M) that could lead to adiabatic conversions. We 
have marked out the projection of the limits that can be achieved by ADMX |30| (Orange) - axion 
search experiment turned into a light shinning through the wall experiment. The bounds put by 
shaded regions with dotted lines come from a summary by [2] and comprise the bounds by both 
theoretical and experimental considerations such as lifetime of the Sun (Red), Horizontal branch 
Star limits (Green), Coulomb law tests (Blue), CMB pollution by the dark photon (Yellow) and 
beam dump experiments E141 and E137 (Purple). In the low mass region the dominant bound 
comes from the drift of fine structure constant (blue, solid/dashed). 
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